FROM SHEAVES ON P 2 TO A GENERALIZATION OF THE 
RADEMACHER EXPANSION 
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Abstract. Moduli spaces of stable coherent sheaves on a surface are of much interest for 
both mathematics and physics. Yoshioka computed generating functions of Poincare poly- 
nomials of such moduli spaces if the surface is P 2 and the rank of the sheaves is 2. Motivated 
by physical arguments, this paper investigates the modular properties of these generating 
functions. It is shown that these functions can be written in terms of the Lerch sum and 
theta function. Based on this, we prove a conjecture by Vafa and Witten, which expresses 
the generating functions of Euler numbers as a mixed mock modular form. Moreover, we de- 
rive an exact formula of Rademacher-type for the Fourier coefficients of this function. This 
formula requires a generalization of the classical Circle Method. This is the first example 
of an exact formula for the Fourier coefficients of mixed mock modular forms, which is of 
independent mathematical interest. 



1. Introduction and Statement of Results 

In the past interactions between physics and mathematics have led to many interesting 
results. Motivated by strong-weak coupling duality (or ^-duality) in physics, this article 
considers various generating functions which appear in the study of moduli spaces of stable 
coherent sheaves on the projective plane P 2 . We express the generating functions of Poincare 
polynomials of moduli spaces of rank 2 sheaves in terms of the Lerch sum and theta function 
which we will recall later. Using these expressions, we prove a conjecture by Vafa and Witten 
[4"5] for the generating functions of Euler numbers. These functions appear to be related to 
Ramanujan's mock theta functions and therefore transform almost as weakly holomorphic 
modular forms, i.e., meromorphic modular forms whose poles (if there are any) may only 
lie in cusps. Our second main result is an exact formula for the Fourier coefficients of these 
generating functions that formally resembles the Rademacher expansion for the coefficients 
of weakly holomorphic modular forms. 

Moduli spaces of coherent sheaves on a complex surface S receive much attention (see for 
example [28J for an extensive work on such moduli spaces). More specifically, one is interested 
in the moduli space Ai(r,Ci,C2) of semi-stable sheaves of rank r with first Chern class c\ 
and second Chern class c%. We will consider topological invariants of Ai, in particular 
the Poincare polynomial p(Ai,s) := ^1=0^ s * an< ^ the Euler number x(A"f) := 

p(Ai,—l), where bi(M) is the zth Betti number: bi(Ai) := dim Hi ( Ai , Z) . Ellingsrud and 
Str0mme [21] computed the Betti numbers of the moduli space of sheaves with rank 1 on P 2 
and other ruled surfaces. Gottsche [23j derived the generating function for p(A4(l, 0, n), s) := 



2000 Mathematics Subject Classification. 11F37, 11P82, 14F05, 14J60, 14D21. 

1 



2 



KATHRIN BRINGMANN AND JAN MANSCHOT 



Yli=obi(-M)s' 1 for rank 1 sheaves on a smooth projective surface, and wrote it as an elegant 
product formula (12. 2p . 

Subsequent work by Yoshioka [49j [50] derived the generating functions of Poincare poly- 
nomials for sheaves of rank 2 on the projective plane P 2 (Eqs. (12. lip and (I2.12p ). Recently, 
the generating functions of the Euler and Betti numbers for rank 3 have also been computed 
[30~1 [35] . Closely related are the computations of the Euler numbers of the moduli spaces of 
vector bundles for rank 2 [29] and for rank 3 [48J. 

These generating functions have also enjoyed much interest in physics, in particular in 
the context of strong-weak coupling duality and instanton moduli spaces. The duality was 
first conjectured by Montonen and Olive as a duality of gauge theory. Their conjecture 
claims that gauge theory with gauge group G and coupling constant g has a dual description 
in terms of the gauge theory with gauge group L G (the Langlands dual group) and coupling 
constant 47r/g. If the theta angle 6 is included in the analysis, then this Z/2 group is enlarged 
to the modular group SL 2 (Z) which acts on the complex parameter t = ^ + ^ l E'H by 
linear fractional transformations. The symmetry of gauge theory under this larger group is 
known as S-duality. 

Vafa and Witten [15] have tested S'-duality for topologically twisted gauge theory with 
M = 4 supersymmetry. They showed that for fixed instanton number the path integral of 
this theory equals the Euler number of a suitable compactification of the instanton moduli 
space, which turns out to be the Gieseker-Maruyama compactification of the moduli space of 
semi-stable sheaves whose Chern classes are determined by the instanton data. S'-duality led 
Vafa and Witten [35] (see in particular Section 3 of [15]) to the conjecture that the generating 
function of the Euler numbers transforms as a (weakly holomorphic) modular form with a 
specific weight and multiplier. These properties were later also understood from the point 
of view of M5-branes, see for example [36] . 

The generating functions ( 12. lip and (I2.12p allow a precise test of the conjectured mod- 
ular properties for rank 2 and P 2 . In fact, we derive similar modular properties for these 
generating functions of Poincare polynomials as for those of Euler numbers. This is quite 
remarkable since present discussions in the literature are limited to the Euler numbers. To 
make the modular properties manifest, we express in Proposition 12. 21 these functions in terms 
of automorphic functions, in particular the Lerch sum (I2.15p . Since the modular properties 
of the Lerch sum are well established, thanks to Zwegers' thesis [51], it is straightforward to 
derive the modular properties of these generating functions. 

Specialization of the Poincare polynomials to the Euler numbers requires one to take the 
derivative of the Lerch sum. Using this relation we prove that the generating function of the 
Euler numbers contains the generating function of the Hurwitz class numbers. A connection 
between the Euler numbers of the moduli space of vector bundles and class numbers was 
earlier proposed by Klyachko [29] . To state our result, let H(n) be the Hurwitz class number, 
i.e., the number of equivalence classes of quadratic forms of discriminant — n, where each 
class C is counted with multiplicity l/Aut(C). We note that H(0) = — ^ and H(3) = |. 
Moreover, we let (throughout q := e 2mT ) 

oo 

(1.1) h 5 {T) := + 3j>" + S j e {0, 1}. 

n=0 



FROM SHEAVES ON P 2 TO A GENERALIZATION OF THE RADEMACHER EXPANSION 



3 



In Section 2, we prove: 

Proposition 1.1. The generating functions of the Euler numbers x (-A4(2, C\, c^j) take the 
form: 



q^Y, X (M(2,-l,n))q n 



oo 
1 v— V 



3fei(r) 
3/i (r) 



t? 6 (t) V(2r) 



q-^ X (M(2,0,n))q n 

n=2 

where r/(r) := Il^Li (1 — <Z n ) ^ s Dedekind's eta-function. 

This proposition is the conjecture [45] mentioned in the abstract. Results of Refs. [29| 14*9] 
led Ref. [45] to this conjecture, and it was verified by a comparison of the first coefficients 
of 3/ii (r)/r/ 6 (r) with Eq. (I2.1ip . The good modular properties of hj(r) after addition of a 
suitable non-holomorphic term (see Eq. (12. 10p ) was a strong confirmation of the S-duality 
conjecture. We refer to Sec. 4.2 of Ref. [45] for more details. 

In the following we recall in more detail what is known about modularity of generating 
functions of class numbers of imaginary quadratic fields. Recall that the Fourier coefficients 
r{n) of 6q, with 0q(t) := Xlnez'?™ > themselves encode class numbers. To be more precise, 
by a result of Gauss, we have that 

(1.2) r(n) = 12 (H(4n) - 2H(n)) . 

If one wants to study the full generating function for the Hurwitz class numbers, then one 
has to move to the world of harmonic weak Maass forms [14J (see also Section [2]). These are 
generalizations of modular forms in that they satisfy the same modular transformation laws 
but instead of being meromorphic they are annihilated by the weight k hyperbolic Laplacian. 
To be more precise, Zagier [51 J showed that the generating function 

h(r) := £ H(n)q n , 

n>0 
n=0,3 (mod 4) 

is a mock modular form with shadow @o(t), notions which we will recall shortly. We note 
that the full generating function for class numbers of real and imaginary quadratic fields 
requires one to consider even more generalized automorphic objects [20J. 

Mock modular forms are related to Ramanujan's so-called mock theta functions, which he 
introduced in his last letter to Hardy (see [43], pp. 127-131) in 17 examples. Ramanujan 
stated that these forms have properties which resemble those of theta functions, but are not 
modular forms. The mock theta functions occur on the one hand in a vast variety of papers 
(see for example [21 [4~J [151 123 HZ] just to mention a few), but were on the other hand not 
well understood for a long time since they lack real modularity properties. The mystery 
surrounding these functions was finally solved by Zwegers in his famous PhD thesis |54| 
in which he related the mock theta functions to harmonic weak Maass forms. Placing the 
mock theta functions into the world of harmonic weak Maass forms has many applications: 
for example the first author and Ono proved an exact formula for the coefficients of one 
of the mock theta functions [TT] and explained how to construct an infinite family of mock 
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theta functions related to Dyson's rank statistic on partitions [F2\- Further applications are 
for example a relation between Hurwitz class numbers and overpartition rank differences 
[9], and a duality relating the coefficients of mock theta functions to coefficients of weakly 
holomorphic modular forms |22| [55] . Part of the difficulty of really understanding the mock 
theta functions was grounded in the fact that these functions have a certain hidden com- 
panion, which Zagier calls the shadow of the mock theta function, and without which the 
mock theta functions are not fully understood. These shadows may be obtained from the 

associated harmonic Maass form by applying the differential operator £2-/0 := 2iy 2 ^ k ^= (with 
k = 1/2 and y := 1m (r)) and turn out to be unary theta functions. Mock modular forms 
are then generalizations of mock theta functions in that the associated shadow does not 
necessarily have to be a unary theta function but may be a general (weakly holomorphic) 
modular form. Mixed mock modular forms are functions which lie in the tensor space of 
mock modular forms and modular forms. 
The functions 

(1-3) /i(r) := -g( = J2^)q n -^, J e {0,1}, 

' ^ ' n=0 

of Proposition 11.11 are examples of such forms. In this paper we prove an exact formula for 
ctj(n), which is the first exact formula for coefficients of mixed mock modular forms. This 
result provides an exact formula for x(.M (2, Ci, C2)), which is clearly of interest for both 
mathematics and physics. The proof requires a generalization of the Hardy-Ramanujan 
Circle Method due to the first author and Mahlburg [10] which may be applied to mixed 
mock modular forms. Let us next place this result in its mathematical context. As usual we 
denote by p(n) the number of partitions of an integer n. Recall that Hardy and Ramanujan 
[25| [26] . in work which gave birth to the Circle Method, derived their famous asymptotic 
formula for the partition function p(n), 

p(n) ■ — -eVV^ ( n ->oo). 

Any/3 

Rademacher [41] then subsequently proved the following exact formula 



x 3 ^ A k (n) ( 7iJ24n - 1\ 
P (n) = 2n(24n - 1)^ £ -fl ■ I J . 

k=i ^ ' 

Here iAx) is the J-Bessel function of order £, and AAn) is the Kloosterman sum 



2 V 12 ^ ~ v ' \\2k 

x (mod 24fc) 
c 2 =-24ra+l (mod 24fc) 



where e(a) := e 2ma and Xi2( x ) '■= (^)- ^ n important tool used to prove the asymptotic 
and exact formulas for p(n) is the fact that 

00 ^ 

n=0 > 
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is a weight —1/2 modular form. Rademacher and Zuckerman |4"2| l52| 153] subsequently 
showed exact formulas for the coefficients of generic weakly holomorphic modular forms of 
negative weight. 

The situation is more complicated if one turns to non-modular objects. Let us mention 
Ramanujan's mock theta functions and in particular 

oo oo n 2 

m = £> )9 " := 1 + £ (i + „).(! + £y-(i + «■)» • 

n=0 n=l ' v 

The problem of obtaining an asymptotic formula for a(n) is greatly complicated by the fact 
that / is not a modular form. Dragonette [19] in her PhD thesis confirmed a conjecture of 
Ramanujan concerning an asymptotic formula for a(n); subsequently this was improved by 
Andrews [Tj. Infinite families of further asymptotic formulas were recently proven by the 
first author [6]. Andrews and Dragonette moreover conjectured an exact formula for a{n) 
which was then proved by the first author and Ono [TT] using the theory of Maass Poincare 
series. The authors of |13| obtained more generally exact formulas for all coefficients of 
mock modular forms of non-positive weight. From the above description it becomes clear 
that currently asymptotic/exact formulas for coefficients of modular or mock modular forms 
are well understood. The situation is totally different for mixed mock modular forms. The 
above mentioned methods cannot be applied as the space of harmonic Maass forms is not 
closed under multiplication. The first such example was considered by the first author and 
Mahlburg [10] and is related to so-called partitions without sequences [3], a partition statistic 
that we do not want to recall for the purpose of this paper. Developing a generalization of 
the Circle Method to involve certain non-modular objects, the authors managed to obtain 
asymptotic expansions for such partitions without sequences. So far this is the only example 
of an asymptotic formula for coefficients of forms in the tensor space. It is of mathematical 
interest to find further such examples. We note that due to the more complicated situation 
the authors of [TU] only obtain an asymptotic and not an exact formula. In this paper we 
derive the first example of such an exact formula. 

Turning back to an exact formula for the coefficients Oj(n) of fj, we require some more 
notation. We let for k G N, g G Z, and u G E 

if g ^ (mod 2k), 
if g = (mod 2k). 



ft a\ j- / \ J sinh 2 (^-^) 

M) fk, g (u):={ } 2 k ™[ 



sinh i — 



2 



Furthermore we define the Kloosterman sums 

Kj y e(n, rn; k) := ipji{h, h', k)e 



— { hn +— 



0<h<k 
(h,k)=l 



where ipje is a multiplier defined in (|3.4|) . and h! is given by the congruence hh' = — 1 
mod k). Finally we let 

2*» == \\ h,9 (I) h (p/(4n-(j + l))(l-« 2 )) i 1 ~ du - 
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Theorem 1.2. The coefficients ctj(n) of fj are given by the following exact formula: 



aj(n) = -1 (4n - U + I))" 1 E «^/ 5 (^^(JTT)) 



fc=l 



*!=i 

oo 



fc=i te{o,i} 

-fc<g<fc 
g=t (mod 2) 

The integrals Zfc )S (n) can be estimated using well-known asymptotic formulas for Bessel 
functions and Proposition 5.1 of |10] . 

Corollary 1.3. The leading asymptotic terms of otj(ri) forn — > oo are: 

aj (n) = (in-i - ^n-i + O (n~ 2 )) e 2 ^. 

We want to make two remarks concerning Theorem 11.21 and Corollary 11.31 

(1) Firstly we note that in contrast to mock modular forms, the shadows of the mixed 
mock modular forms do contribute to our leading asymptotic terms. One could 
determine further polynomial lower order main terms. 

(2) Secondly, we like to mention that the first term in the exact formula are the coef- 
ficients of a negative weight Poincare series as described by Niebur [40| . Numerical 
experiments by F. Stromberg give strong evidence that generically this term does not 
converge to an integer. G. W. Moore and the second author |33| considered negative 
weight Poincare series (which are essentially a sum over ro^SL^Z)), because of their 
interpretation in the context of the correspondence between 3-dimensional Anti-de 
Sitter space and 2-dimensional conformal field theory [18J. This led them to consider 
alternate functions, say fj, in addition to fj, whose coefficients are given by the first 
term of Theorem 11.21 Interestingly, the theorem shows that fj appears naturally as 
part of fj. Moreover, it is possible to show that fj — fj can also be written as a sum 
over r oo \SL 2 (Z). We leave a precise discussion for the future. 

The results of this paper might be relevant for various applications and current developments, 
some of which we want to list here: 

(1) The appearance of Lerch sums in the generating functions of Poincare polynomials 
is rather intriguing. Besides for P 2 , one can show, using the results of [49, 50J, that 
they also appear for rank 2 sheaves on ruled surfaces. Their appearance is essentially 
a consequence of the contributions of stable bundles to the generating series for a 
specific polarization. It would be interesting to investigate whether Lerch sums play 
also a role for different surfaces, higher rank sheaves, and related systems like Calabi- 
Yau black holes. 
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(2) The exact formula for the coefficients of fj can be generalized to other mixed mock 
modular forms. Besides the intrinsic mathematical interest, this might also prove to 
be very useful in physics, in particular in discussions on black hole entropy and the 
AdS3/CFT 2 correspondence. For example the functions hj/r) 2A are known to appear 
as generating functions of the degeneracies of M = 4 dyons [16J. 

The outline of this article is as follows. Section [2] reviews briefly the generating functions 
of invariants of moduli spaces of stable sheaves, relates those of rank 2 to the Lerch sum and 
theta function, and proves the conjecture by Vafa and Witten. Section [3] derives the exact 
formula for the Fourier coefficients of fj, using the Hardy- Ramanuj an Circle Method. 
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2. Generating functions of topological invariants 

Let us start by recalling some of the relevant background. We refer the reader who is 
unfamiliar with these notions from algebraic geometry to consult textbooks like [241 [39] . It 
is well-known that the Betti numbers of P 2 equal bo = 6 2 = &4 = 1 and b\ = 63 = 0, therefore 
%(P 2 ) = 3. The total Chern class c(P 2 ) of the tangent bundle of P 2 is defined as: 

C (P 2 ) 1=1+0! (P 2 )+C 2 (P 2 ) = (l + J) 3 , 

where J is the hyperplane class and q(P 2 ) G H 2i (F 2 ,7j). Since these cohomology groups 
are 1-dimensional, we will also denote the integrated forms J Cj(P 2 ) by q(P 2 ), thus Ci(P 2 ) = 
c 2 (P 2 ) = 3. 

Chern classes Ci(E) are defined for any sheaf E on P 2 , and play a central role in the 
classification of sheaves. If no confusion can arise, the Chern classes Ci(E) are in the following 
abbreviated by Cj. The complex dimension of the moduli space of stable sheaves may be 
written in terms of these Chern classes as 

(2.1) dim c (M(r, c u c 2 )) = 2rc 2 - (r - l)c 2 - r 2 + 1, 

where r is the rank of the sheaf. The moduli space of stable sheaves generically depends on 
the choice of an ample line bundle over the surface. However, since 6 2 (P 2 ) = 1, stability does 
not depend on this choice. 

We are interested in the generating functions of the Poincare polynomials and Euler 
numbers of moduli spaces Ai(r, Ci,c 2 ) as functions of c\ and c 2 . Twisting a sheaf by a 
line bundle E ® 0(k) gives an isomorphism between the moduli spaces Ai(r, ci,c 2 ) and 
M.(c.\ + rk, c 2 + (r — l)kc\ + |r(r — l)k 2 ). It is therefore sufficient to only consider c\ 
(mod r). 
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The generating function of the Poincare polynomials of .M(1,0, C2) for any surface S is 
given by 
(2.2) 



'Y _|_ s 2m-l^mj 1, i( s ) _|_ s 2m+l^m^ b i( s ) 



m i (1 - S 2(m-l) t m)6o(S) ^ _ s 2m t m)MS) ^ _ s 2(m+l) t m)6 (S) ' 

To exhibit the modular properties for S* = P 2 , we write it in terms of the Jacobi theta 

function 9\(z; r), which has the following sum and product expansion (w = e 2ntz ): 

(2.3) 

2 00 
9x(z; t) := i (-l) r ~^ ¥ ^ r = *?* M - II ( X ~ ( X ~ i 1 ~ w ~ l( l n ) ■ 

reZ+| n=l 

This theta function transforms under the generators S :— (1 ~q) and T := (J }) of SL 2 (Z) 
as: 

= -^A/^irexp ^ ^i(^; T ), 

Moreover #1 has simple zeros at the points z = nr + m with n, m G Z. 

With the substitutions g = s 2 t = exp(27rir) and w = s 2 = exp(2iriz), equation ( 12.2D 
becomes 



(2.4) (-M(l,0,n),w^ (gw -1 )' 



1 

I I «J2 — W 2 



n>0 



0i(«;r) 



The Betti numbers can be obtained by first expanding equation ( 12.4ft in g for g « 0, and 
then in w for w ~ 0. One easily sees that ( 12.4ft has no poles for zGZ, but does have simple 
poles for 2; = mr + n, with (m,n) G Z 2 , m 7^ 0. The Fourier coefficients of ( I2.4[) depend 
therefore on the choice of contour to extract the Fourier coefficients. The physical origin of 
these poles is however unclear, but might be related to the fact that the Poincare polynomial 
is not a supersymmetric index like the Euler number. 

The above substitutions for s and t are not arbitrary but compatible with the Lefshetz 
s/(2)-action on the moduli space. If J3 is identified with the Cartan element sl(2), then the 
action of J 3 on an harmonic form on the moduli space is given by |17| |2"4"] 

J 3 cu = ~ (deg co — dim M) uj. 

The eigenvalue of J3 is the exponent of w in the expansion. 

The Euler characteristics are obtained by setting s = —1: p(Ai(l, 0, n), —1) = x(.M (1, 0, n)). 
Then the generating function becomes: 



00 

(2.5) / 1>0 (r) := q-^x{M(l t f n)) q n = 



n=0 
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thus a modular form of weight —3/2. For r > 1, one can more generally define the functions 
(2.6) fr, ci (r):= £ X (M(r, c 1; c 2 )) ^a-^)/^ 



,r— 1 2 
■ 2r C l 



with A the discriminant of A := - (c 2 — ^c 2 ) and < Ci < r — 1. These functions are 
expected to exhibit transformation properties of a vector-valued modular form of length r 
and weight — x(P 2 )/2. The modular properties of this function are most straightforwardly 
derived from the point of view of multiple M5-branes wrapping P 2 £g> T 2 . This leads to a 
generating function which also sums over all C\ |36[ 



(2.7) Z r (p;r):= J2 x(M(r,c u c 2 ))q r y A ~' 24 'J q 2r V i: 2 J(-^+^t^ 

Cl ,C2&Z 

with £ := e 2mp . Physical arguments suggest that this function transforms under SL 2 (Z) 
like a Jacobi form which is non-holomorphic in r and has weight (§, — §)■ Moreover, the 
isomorphism of moduli spaces due to twisting by a line bundle implies a decomposition of 
Z r (p\ r) into theta functions O r . M (p; r) and vector- valued modular forms / r , M (r): 

(2-8) Z T {p-r)= V / rjt (r)e v (p;r), 



(mod r) 



with 



8 r >;r):= ^ g 2r V 2 ^ (-£)' 

n=fi (mod r) 

This decomposition implies that 

D r {Z r {p; r)) = 

with D r := J-p + 4^7 ^p-. Functions which satisfy this condition together with an appropriate 
transformation law are known as skew (weakly) holomorphic Jacobi forms (Hj. In particular 
for r = 1 we have that 

Zl{p]T) = . 

Already for rank 2, we will find two refinements of these physical expectations. As ex- 
plained in the introduction, the /^(t) appear to be mixed mock modular forms, such that 

D 2 (Z 2 (p;r)) + 0. This is in physics called a "holomorphic anomaly". The other refinement 

concerns the integrality of the Fourier coefficients. 

Before returning to the functions of interest for this paper, we want to recall the precise 
definition of harmonic weak Maass forms. Here we only require the case of half-integral 
weight on T (4). 

Definition 2.1. A harmonic weak Maass form of weight k G | + Z for the group To(4) is a 
smooth function / : H — > C satisfying the following: 
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(1) For all ( a c b d ) E r (4), we have that 



Here (|) denotes the Jacobi symbol, = 1 for d = 1 (mod 4) and = i for d = 3 
(mod 4), and yfr is the principal branch of the holomorphic square root. 

(2) We have A^f = 0, where (r = x + iy) the weight k hyperbolic Laplacian is defined 

as 

Afc ■= -y 2 (— + —) + ik (— + I — 

\dx 2 dy 2 J \dx dy 

(3) The function / has at most linear exponential growth at all the cusps. 
Using this notation, the function 

;2. 9 ) Mr) := £ W + il±il f -^W d „ 



16tt 

n = 
n = 0,3 (mod 4) 



T + W 2 



is a harmonic Maass form of weight § on T (4) (see [51J). We moreover require the restrictions 

of h to arithmetic progressions 0,3 (mod 4) (individually). It is not hard to see that the 
associated harmonic weak Maass forms are given by 

(2.10) h (T):=h {r)+ (l±A n^i^dw, JG{0,1}, 

8vr J_- ( r + U7 )2 

where the functions hj were defined in ( 11.11) and Qj(r) := gi( 2ra+ j) . 

nez 

Now we continue our discussion on generating functions related to semi-stable coherent 
sheaves of rank 2 on P 2 . Yoshioka [4~9|[50] computed the generating functions of the Poincare 
polynomials. To present his result, define 

Zs := (l-t)(l-st)(l-sH)' 

The expression for the generating function for c\ = — 1 is [49J: 

~ n rf>1 z s2 (F 2 ,s M -H d ) 2 



n= l V I L^n^L 

s 2( 6+l)(26+l) g 2b(2fr + 5) ^ . 

3 8(6+l)£26+l ~~ l _ s Sbf2b 



b>0 



x j2 ( - J ^ .... - •" . .. . ) 1 
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and similarly for c\ = |50| : 

(2.12) J2p(M(2, 0, n), s) f> = (1 _ 



n=2 

„2(6+l)(2&+3) „26(26+7) \ „2(6+l)(26+l) _ „26(26+l) 



/ „2 (6+1 26+3 „26 26+7 \ 

/ s ^ s \ 6 2 +3fc+1 

\1 — s 8(6+l)^-26+l I _ s 86^26+l ) ~ r 2s 2 

v 6>0 V 7 6>0 



t 6(6+l) 



2s 2 (1 + s 2 ) 

Here we have corrected a sign error in Remark 4.6 of |50| . 

To simplify the expressions ( 12. lip and ( 12. 12ft . we make the substitutions sH = q and 
s 2 = w, analogous to the substitutions in the rank 1 case. One finds after a straightforward 
computation: 

Proposition 2.2. The generating functions of the Poincare polynomials p (Ai(2, C\, C2), s) 
take the form: 



00 

(2.13) q-s^p (M(2,-l,n),w^ (qw~ 2 ) n 



n=l 



(1 W) ti(2z-r,l-T-z;2r) 



W2 Q\{z\ t) 
00 

(2.14) q~*^2p (M(2,0,n),w^ (g^ 2 )" = 

n=2 

(l-w)fl _i 3 / 1 N \ iy) 

n „„. o u 2^ — r, 7T — 2; 2r — — -. 

w 2 6 2 (z;t) V2 V 2 V 2w 2 ^(2z;2r) 

where h(u,v;t) is the Lerch sum defined by 

^ttiu ^ / ^\n^,ni(n 2 +n)T+2irinv 

It, D G C. 

Moreover, we define 

(1 — w) q~^ q n2 w~ n 



f (7-t\ = [ ' q 4 V 

/2 ' U ' J • ^ 2 (^r)^;2r)^l-g 2 ™-% 2 ' 



n£2 

3 5 



(1 — iu) / 1 q 4.W2 s—^ q n +n w n 
hfi ^ T > := w i e 2 {z; r) \ 2 + 6 2 {z; 2r) ^ 1 - q^w 2 

with 
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and 



91(Z;T) '= 9^) ^ 1 - g 2n-i w2 = -A* (2^ - r, 5 - r - z; 2r) , 

1 q'iwi v-^ qn 2 +n w -n 1 _1 8 / i -. 

: — I - tt^ r > n — ; — o = Q i w 2 u (2z — t, ^ — z;2r) 

2 # 2 (z;2r) ^ l-g 2 "" 1 ™ 2 2 H PV ' 2 ' ; 



Similarly to the case of r = 1 these functions have poles for z G mr + n with (to, n) G 
Z 2 , m ^ 0. 

Since we have now explicit expressions for the generating functions of Poincare polyno- 
mials for r = 1,2 at our disposal, it is particularly interesting to investigate their modular 
properties. For rank 1, equations (12.41) and (12. 5p show that the generating function of Euler 
numbers is indeed a weakly holomorphic modular form, whereas the generating function for 
Poincare polynomials transforms as a Jacobi form of weight — | and index — | (up to the 

pref actor — w~z). 

To make the modular properties of ( I2.13P and ( I2.14p more manifest, we recall some results 
of Zwegers' thesis |54j . The Lerch sum ( I2.15P does not transform as a Jacobi form under 
SL 2 (Z). However, the completed function 



fi(u, v; t) := jj,{u, v; t) + 2^( M — v 'i T ) 

transforms as a mult i- variable Jacobi form of weight |. Here the function R(u\ r) is defined 
by: 

r— = „—irir 2 r—2niru 



R(u;t):= (sgn(r) -E ((r + a)\/2y)) (-l) r ~V 

T 

with a := lm(u)/y, and 



E(z) := 2 / e^du. 
Jo 

To be more precise, we have that: 

(1) For k, l,m,n G Z, we have that: 

ji(u + kT + l,v + mr + n;r) = e m(*-m) a T+MMM^ tti V] r ). 

(2) For 7 = (» J}) G SL 2 (Z), we have that: 
u f ar + &\ , 3 n i _7ricCu-^ 



// — , — ; — )=u(7) (cr + rf)2 e fi(u,v;r) 

ct + a ct + a ct + a ' 



with «( 7 ) := r? (Sg) / ((ct- + rf)^(r)) . 



Moreover we require the following identity, which allows us to shift parameters in the function 

H (z G C): 

2.16 /i (u + z,v + ^;r - /i (u,v;t = — — — — ■ -. 

0i(u; r)0i(u; r)0i(ii + z; r)6>i(i; + 2; r) 
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We now turn back to the functions gj(z; r) and define their completions as: 



gj(z;r) := &(z;t) + -Rj(z;r) 



with 



We note that 



Rj{z; t) := ( s S n ( n ) " E ^ 2n + 3a ) V^)) Q 



R (z;t) = — 1 — iq~*w%R ^3z — r — ^; 2r 

Using the above stated transformation properties of ju, one can show that the functions ^ 
are invariant under T 4 , and transform under S~ 1 T~ 4: S G To(4) as 

* = (4- + 1)* exp (2vr^) r). 

Moreover, one can prove with some more work that the function g<j(z\ r) may be viewed as 
components of a function that transforms like a vector valued modular form for SL 2 (Z) of 
weight | and with the same multipliers as the hj(r) defined in equation (12. 10p . 

We observe that if f2,j(z; t) would be completed to f2j(z', r) by changing gj(z; r) to "gj(z; r), 
they would transform as Jacobi forms for r (4) of weight — \ and index —^f, if we ignore the 

prefactors (1 — w)/w 2 ~ J '. The non-holomorphic parts of f2,j(z; r) might appear naturally in 
physics, but precisely how is unknown. Since the functions gj{z; r) transform as a modular 

vector, the function Z 2 (z,p;r) = ^j=o i f2,j(z\ t)^2j(p\ t) transforms with weight (|, — |) 
under SL 2 (Z) (ignoring the prefactors), which can be understood from physics. 

As expected from physical arguments, the modular properties improve if one takes the 
limit »2 — > — 1. One can derive straightforwardly that in this case the last term in equation 
( I2.14p is equal to ir] _3 (2r), which is a modular form of To (4) with a non-trivial multiplier. 
More interesting is that the limit translates to taking the derivative of the Lerch sums in Eqs. 
fl2~T3|) and (l2TT4j) . Proposition O gives for these terms f 2 j(r) : = / 2 j(0;t) = 3hj{r)/rj 6 (r). 

Before proving the Proposition 11.14 we would like to make a couple of remarks concerning 
/2j(t), and 2 2 (p; t) defined by equation (|2.7|) . The completions fzj(T) can be obtained from 
f2,j(z; r), by computing the coefficient of z 1 in the Taylor expansion of cjj(z; r). Due to the 
non-holomorphic term, D 2 (^Z 2 (p; r)J 7^ 0. One finds 

^ r 2(P ' T) J = 16^/2 ^ (r) ' 

which is proportional to Zf(p;r). The authors of |36| conjecture that such an anomaly 
appears generically for r > 2. We did not find such a factorization in the case of Poincare 

polynomials, that is to say for D2 [Z2(z,p; r) 
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The generating functions of the Euler numbers, /2,i( r ) and /2,o( r ) + \v~ 3 {^ T )^ do not com- 
bine to a vector- valued modular form because of the term |r/ _3 (2r). This term disappears, 
if we consider the generating functions of the rational invariants 

X(r):= (-l) dimc(A " (r/m)) x(r/m)/m 2 , 

m>l,m\T 

where T represents the data of the sheaf (r, c±, C2). Refs. [34, 45J give also evidence that the 
generating function of the rational invariants have better modular properties, than the 
ones for the integer invariants This is the second refinement, alluded to below equation 

On the other hand, the coefficients of /2,o( r ) + jV~ 3 (^ T ) are required to be integers. This 
can easily be seen from the arithmetic properties of the functions. To see this, multiply the 
function by ^(t) 6 , which gives 3/io(t) + |Oq(t + \). Integrality of the coefficients of this 
function is manifest, due to the properties of the class numbers H{n) and 0q(t + |). 



Proof of Proposition 11.11 One could prove the proposition straightforwardly by verifying 
1) that the shadows of f2 s j{j)i as obtained from f2,j(z', t), coincide with those of 3/ij(r)/ry 6 (T), 
and 2) that a specific number (related to the dimension of the space of associated modular 
forms) of coefficients agree. Since this is rather technical, we choose to prove the proposition 
by relating it to known expressions in the literature. 

We start with the identity for / 2 ,i- The limit z — > of f2,i{z',r) is finite and leads to 
differentation of the Lerch sum: 



(2.17) 



d 



?7 6 (r) dw 



H I 2z — r, 



t+- 2 ,2t 



w=l 



Using fT2TT5]) yields that 



fl I 2,z — T, 



t+- 2 ,2t 



1 1 

+ 



-t,-3z-t+ -;2r 



zV (2r)0i (-z - 2r + §; 2r) e x {2 Z ] 2r) 



9 1 (2z-r;2r)9 1 (- 



Z — T 



|; 2r) Q x (-T- 2r)e x (-32 - r + §; 2r) ' 



One can prove that the second summand contributes \®\{t) / rf {r) to f )2.17p . Moreover, the 
contribution from the first summand is given by 



Q 



Using work of Kronecker 
(2.18) hxir) 



2 V 6 (r) o (r 
Mordell 
1 



^(2n-l)g 



2n-l 



and Watson 



2O (r) 



Q 



i^ (2n-l)g i 3 



one can prove that 
1 
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From this the claim may be easily concluded. We first note that by Watson (correcting a 
typo) we obtain that 



(2.19) 



n=0 6K ' neZ 



qi n — q n 2 



where F(n) counts the number of uneven equivalence classes of positive definite quadratic 
forms of discriminant —n. Next one can easily show (for example by using the theory of 
modular forms) that 



ef(r) = ^r(4n + 3)<f 

n=0 

where the coefficient r(n) is defined by 

oo 

e*(r) = 5>(n)<f . 



n=0 



Now a direct computation gives (I2.18p . 

We next turn to ^q. It is not hard to see that 



(2.20) / 2 , (r) 
We find using (j2"TT6"|) 



1 d 



rf{r) dw 



q *w 2 fi [ 2z — r, — z + -; 2r 



- w=l 



li[2z-T,-z+ -;2r 



-3z + r; 2r 



ir ) 3 (2r)e 1 (-z; 2t)0i (2z - r + ±; 2r) 



d^z - r; 2r)0i (-z + §; 2r) X (-|; 2r) <?i(-3z + r; 2 



T 



One can show that the contribution of the second summand to (I2.20p equals — / rf (t) . 

Moreover one can prove that the first summand gives a contribution of 

n(-l) n q n2 



V 6 (r) O (r + ±) 



+ q 



2n 



Now the claim easily follows using ( II. 2p and the identity 



(2.21) 



5> 

n=0 



V + ^l 2 - i 



n(-l) n <f' 



20 o(r+|)^ 1 + ^ 12 ° 



e 3 n (T). 



Indeed, equation (I2.2ip may for example be concluded by combining Theorem 1.1 and Corol- 
lary 1.6 of [8] and inserting the generating function for / given in [9]. □ 
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3. Exact formulas for a,-(n) 

The introduction motivates the derivation of an exact formula of the Fourier coefficients 
of /o( r ) = |/2,o(t) and /i(r) = |/2,i(r). This will be the subject of this section. We start by 
providing various useful transformation formulas, after which we use the Hardy- Ramanuj an 
Circle Method, to derive the exact formula. 



3.1. Some transformation formulas. In this section, we give transformation properties 
for the class number generating functions h and hi. Throughout, we let z G C with 
Re(z) > 0, k > 0, (h, k) = 1, and h' defined via the congruence hh' = —1 (mod k). 
Moreover, we assume that A\h' if k is odd. We require the transformation law of the eta- 
function: 



(3.1) 

where uj h ^ k is given by 



rj [ -(h + iz) 



;{h-h') 



-1 
J h,k 



_ 1 
Z 2 



-'/,./, := cxp j 7li ^2 

(mod k) 




with 



, x — [x\ — \ if x G R \ Z, 
^ X)) 1 if x G Z. 



Writing O and ©i as eta-quotients 



©o(r) 



7? 5 (27 



1 (r) = 2 



t? 2 (4t) 



r7 2 (r)r7 2 (4r)' ~ ±vy ~ r](2r) 

yields the following transformation law (j G {0, 1}): 

• ^ =4 E xAh,h',k)Q e (Uh'+ t - 
v te{o,i} V V 

Here the multipliers Xjt are defined as follows: 



Xoo(^, h',k) :-- 



Xio(h,ti,k) :-- 



,2 , ,2 



V2 





, fe 

, ,2 
' J h,k UJ 4h,k 

~77> 

u 2h,k 







e 2*; 



1 ^2h,k 
^-e 2fc 



2ft 



if 4|fc, 

if 2 f fc, 
if 2||fc, 

if A\k, 
if 2 f k, 

i£2\\k, 



i 



X oi(h,h',k) := { V2 



,2 , ,2 
u h,k w Ah.k 



2 2 

^h,k^ 2h k nih > 
-e 2k 



^7 



Xu(h,h',k) : = 



7i 



4h,fc 



if A\k, 
if 2 f k, 

i£2\\k, 

if 4|fc, 

if 2 f Jfe, 
if 2||fc. 
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We next use the well-known behavior of hj under inversion and translation 



17 



h {r + 1) = ho(r) h 1 (r+ 1) = -ih x (j) 

*• ( ~) = (Mr)+Ur)) *, = (Mr) - S.M 

This gives that /i,- has multiplier dual to the one of 0,. To be more precise we have 



hj[\(h + iz)\=-z i Xjt(h,h',k)h t (^(h' + 
' te{o,i} ^ ^ 



From this, a straightforward calculation shows that 



(3.2) ( ^{h + iz)j =-z~ 3 z Xn(h,h',k)h e (^(h'+- 
' te{o,i} ^ ^ 



4V27T Jo {t+ ±) 

Defining 



[W + X)2 



Ij{x) := / — -F-tfou 



we may rewrite ( 13.21) as 

h j (-(h + iz)\=-z- 3 z Xjt{h,h',k)ht(^- + 

Dividing by 77 s and applying ( 13.1 ft yields that 

(3.3) + =.te^-^ + 



te{o,i} 

h' TTt(j + l)fe 



te{0,l} 
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with 
(3.4) 



i/j 00 (h, h',k) 



1 u h,k^2h,k 
V2 "\ 



i/)ia(h, ti, k) := < 



V? ^2h,k 



if i\k, 

if 2 f k, 
if 2 1| /fc, 

if 4|Jb, 
if 2 f k, 

if 2||Jfe, 



ip 01 (h,ti,k) :-- 



1 u >h,k^2h,k 
V2 



if 4|fc, 
if 2 f fc, 



i — 2fc if 2 fc, 



(/i, ti,k) 



-s-^e 2* 

1 w h,fc w 4h,fe 







if 4|Jfe, 

if 2 f k, 
H2\\k. 



For later purposes we require a different representation of Ie(x). Similarly as in [8], one can 
show 



Lemma 3.1. We have for igC Re(x) > 



(3.5) 



E 



g (mod 2fc) 
g=j (mod 2) 



4fc 



y/x y/2irk 2 X j - x 



-27TIEM 2 



where Sq^ = unless g = (mod 2fc) m which case it equals 1. 



Note that we corrected a sign error in the statement of Lemma 4.4 of j8]. 

3.2. Proof of Theorem II. 2L Throughout this section, we use the notation from Subsection 
13.11 For the proof of Theorem 11.21 we employ the Hardy- Ramanuj an Circle Method [42j and 
write for j G {0, 1} 

oo 

j+1 



n=0 



fM) := g^/i(r) = 5>(n)<f . 
By Cauchy's Theorem we have for n > 

otj{n) = 



, 4k 

2vri J c q n+l 

where C is an arbitrary path inside the unit circle looping around counterclockwise. We 

-2tt 

choose the circle with radius r = e"^, where we later let iV — > oo, and decompose it into 
consecutive Farey arcs of order N: 



aj(n) 



E 

0<h<k<N 
(h,k)=l 



f I ^-+2iri£+2ni( l 
fj * 



e 



-2nin<f> 



with: 
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where ^ < | < |j| are adjacent Farey fractions in the Farey sequence of order N. From the 
theory of Farey fractions it is known that 

' <tAt (7 = 1.2). 



k + kj ~ N +1 

Using the transformation law ( 13. 3 h and z = k(^ ~ ^0), we obtain 

// 

i"9 / 1 / 

■Kih! h,k I X / 2 \ \ 2irz/„ j + 3 



e^^~ 2^ ip je {h,h' ,k)e^ ^ M fc I 

0<h<fc<]V £ g ( 1} ■'-' d h.k ^ ^ 

(h,fc) = l L J 

4V27T n ^rr^« ^-1 



(fc,fc)=l 



e k \ n " 4. ) Z 2 



We will abbreviate the first summand by ^ 1 and the second by J^ 2 . 

We first consider ^ 1 and split of the terms with negative exponent in the Fourier expansion 
as they contribute to the main term. For this, we write 



fo{r) = -—q * + ^&o(n)g n i 

n>0 

/i(r) = J>(n)<T*. 



n>0 



We denote the contributions of the negative exponent to by Using the estimates 

v'h.kiK.k < m and N 2 < N ~ 2 § ives that 



1 1 \ 0<h<k<N / 1 



To estimate J^ 2 , we use the representation of Xg given in Lemma 13.11 We start with the 
contribution of the first summand in the representation of Xg (only occuring for i = and 
g = (mod 2k)). Splitting of the non-principal terms yields as before an error of N~2. To 
estimate the remaining terms of J^ 2 , we aim to estimate integrals of the shape 



/°° 2 2 

e~^fk, g (u) du. 
oo 



Similarly to the case of Fourier expansions we are interested in the "principal integral part" 
contribution. To be more precise, we let for b > and g G Z, 



2-iri) 5 



"V^ „ 2 
^ lj _ . 2ttu 

Jk,g,b{z) := e fc* I e kz f k>g (u) du. 
Vb 



Similarly as in [10J, we may show: 

Lemma 3.2. As z — > oo we have for —k < g < k: 
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(1) lfb<0, then 
(3.6) 



\Zk,g,b{z)\ < W X { Y 



11 if 9 :' C 
1 ifg = 0. 



(2) lfb>0, then 

Zk,g,b{z) = Jk,g,b{z) + £k,g,b{z), 

where the error £k,g,b satisfies the same estimate as Ik,g,b i n Ii3.6\) . 
Proof. Recall that Re Q) > §. We use the estimate 



sinh 



ixu mg 
~k~^k 



cosh 



TXU 

~~k 



— 7TI 



g_ 1 

2k 2 



> 



sin 



T^g_ 

2k 



Now for < x < ^, the function j s bounded from below, thus for —k < g < k, g ^ 



Moreover 



We now define h g ^ as 



1 



sinh fx) x 2 



1 



1 



x 2 sinh (a;) 



< 1. 



h 



,,2 



g,k 



if -k<g<k,g^0, 



Then by the above 
We now first assume that b < 0. Then 

POO 

\Ik,g,b(z)\ < \ z \^ h g,k / e 



1 if 5 = 0. 

\fk,g( u )\ < ^g.fc- 



dit < |z| 2 hg,k\ TTTT\ < M 2 ^?,fc 



Mi) 



which gives the claim for b < 0. The case 6 > works similarly. 



□ 



The contribution of non-principal part of the remaining terms of may be estimated 
against a constant times 



1 r^h h 

J2k / Sl X ^o(^)|rf0« 



/(..fc 



y: 1 



h,k g 



h,k 



k J-d 



\Z\ 2 



"5? 



In the terms coming from the principal part, we may similarly truncate the integral to lead 
J kg i(z). Combining the above, we have shown that 



a j (n) = S 1 + S 2 + S 3 + o(N- 
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with 



i rfi" 

1 ^— \ -2irihn , / M ggz/ _ (j + m + _n_ 2 

51 := - — > e * ip jQ (h,h',k) e k ^ * J+ 2feZ2;2 

0<h<k<N 17 "fc.fe 

(h,fc) = l 

1 ^ — -\ r— -2irihn P h,k gjr 

5 2 := —7=- > Vfce * if> jQ (h,h',k) e« 



2 I z(' j + , _ti_ 2 



0<h<k<N 
(h,k) = l 



-6 £ 

Q<h<k<N ^6{0,1} V y h.fe 



8vr 

1<A<A<JV te{o,i} 
(/l ' fe)=1 -k<g<k 

g=£ (mod 2) 

We next write the path of integration in a symmetrized way 



i 



h.fc kN kN k(k+k 2 ) 

The second and third term contribute to the error term and may be estimated as before. To 
finish the proof, we require estimates for integrals of the form (r > 0) 



i 

" kN 

-Ek,r,n,m • / Z C 
"kN 



In a standard way (we refer the reader to [5] for the details) one may show that 

2tt p\f T (Aix . \ ^ ( 1 

-Lk,r,n,m = ~r ~ J *r+l -rV nm + O 



k \nJ V k J ykN^ 1 

Inserting this bound into the Si and letting N — > oo now easily gives the claim. □ 
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